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In a recent paper, the generalization of the Jensen Shannon divergence (JSD) in the context 
of quantum theory has been studied (Phys. Rev. A 72, 052310 (2005)). This distance between 
quantum states has shown to verify several of the properties required for a good distinguishability 
measure. Here we investigate the metric character of this distance. More precisely we show, formally 
for pure states and by means of a numerical procedure for mixed states, that its square root verifies 
the triangle inequality. 
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I. INTRODUCTION 



Fundamental physical theories are formulated in terms of an abstract space. This is the case of Relativity Theory, 
Quantum Mechanics (QM), Yang-Mills like theories, and every proposal for Unified Field Theory. On each abstract 
space different structures can be defined. For example topological, differentiable, affine and metric structures are 
ubiquitous in space-time models. A prescription for measuring just how close two points of the concomitant space 
are is what we mean here by a metric structure. A more precise distinction between a distance and a metric will be 
given below. 

In principle each one of above mentioned structures can be defined in an independent way. In Special Relativity 
theory the space-time is the standard manifold provided with the (fixed, non-dynamical) Minkowskian metric. In 
General Relativity, instead, the space-time is a differentiable four dimensional manifold where the metric is given by 
the matter-energy distribution (throughout the Einstein's field equations). In both cases the metric is compatible 
with Lorentz's covariance. It is worth mentioning here (as known since the pioneering works of Gauss and Riemann) 
that the metric defines every geometrical property of a differentiable manifold. 

In QM the corresponding abstract space is a (finite or infinite dimensional) Hilbert space Ti.. In its mathematical 
formalism the states of a physical system S are represented by operators (density operators) acting on Ti.. More 
precisely the states of the system S are represented by the elements of B{H)i , that is, the set of positive trace one 
operators on H. The notion of a state as a unit vector of H refers to the extremal elements of B{Tl)i {p e B{H)f is 
extremal if and only if it is idempotent, = p). In this case p is of the form for some unit vector \ip) sTi, and 

is called a pure state. 

In the case of a Hilbert space, the basic underlying structure is that of a vectorial space provided with an internal 
product ( I ) between elements of Ti. From this inner product several ways of measuring "proximity" between two 
elements of Ti, can be defined. For example, the Wootters's distance 

dw{\v), m = dw{W){vl IV-XV-I) = arccos(|(^|^)|) (1) 

is a very important one. On one side ([T]) represents the angle between the (pure) states \lp) and \il))] on the other, it 
has to do with the statistical fluctuations in the outcomes of measurements into the QM formalism [l| . Finally, ^ 
is invariant under unitary evolution. Therefore, we can think of ([1]) as a very natural distance between pure states in 
QM, in some sense imposed by the quantum theory itself. A generalization of this distance to mixed states have been 
studied by Braunstein and Gaves 

Before going on let us remind the reader of a formal distance-definition. Let X be an abstract set. A function 

d : X X X ^ R 

is a distance defined over the set X, if for every x, y eX, it satisfies the following properties: 



d{x,y) > for x^y and d{x,x)=0 
d{x,y) = d{y,x) 



(2) 
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If, for every x, z eX, the function d also verifies the triangle inequality: 

c?(x, y) + d(y^ z) — d{x, z) > (3) 

it is said that c? is a metric for the space X. Incidentally, we mention that the function given by ([!]) is a metric. 
However, only a few among all distances between quantum states historically introduced in the literature verify 
condition ([3]). 

The definition of distance between mixed quantum states is a topic of permanent interest. This interest has been 
lately rekindled on account of problems emerging in information theory (QIT) i, H i, i, 13 ■ 

In introducing distances 

between quantum states, different roads have been traversed. We have already mentioned the case of the Wootters's 
distance and its generalization, presented in [^. Recently, a rather interesting approach has been advanced by Lee 
et al. in reference Q. There these authors characterize the degree of closeness of two states with regards to the 
information that can be attained for each of them from a complete set of mutually complementary measurements plus 
an invariance criterium. The resulting distance-measure is equivalent to the Hilbert-Schmidt metric. Let us recall 
that this metric emerges from the primitive structure of the Hilbert space. Indeed, an inner product between bounded 
operators acting over the Hilbert space H can be defined in the fashion 

{A\B) = Tr{A^B) 

The Hilbert-Schmidt norm of the operator A is given by — {A\A) and from this, the Hilbert-Schmidt metric 

between two operators A and B is defined as 

dHs{A,B) = \\A-B\\Hs (4) 

Another way of dealing with the problem of introducing distances between quantum states is generalizing the 
notions of distance defined in the space of classical probability distributions. This is the case of the relative entropy, 
which is a generalization of information theoretic KuUback-Leibler divergence. The relative entropy of an operator p 
with respect to an operator cr, both belonging to B{Ti)i , is 

S(p, a) ^Tr[p{\ogp- Yoga)], (5) 

where log stands for logarithm in base two. The relative entropy is not a distance (and obviously is not a metric 
either) because it is not symmetric and does not verify the triangle inequality ([3]). Worst, it may even be unbounded. 
In particular, the relative entropy is well defined only when the support of a is equal to or larger than that of p Q 
(the support of an operator is the subspace spanned by the eigenvectors of the operator with nonzero eigenvalues). 
This is a strong restriction which is violated in some physically relevant situations, as for example when cr is a pure 
reference state. 

To overcome such problems we have recently investigated an alternative to the relative entropy 8| that emerges as 
a natural extension of a symmetrized version of the KuUback-Leibler divergence to the realm of quantum theory. In 
the classical context this quantity is known as the Jensen-Shannon divergence (JSD) and was introduced by C. Rao 
^and, indegendently, by J. Lin [l^]- It has been applied to a diversity of problems arising in statistics and physics 
11 , [l^, [l^ O, [l5 L Among its most significant properties one can include its boundedness and its metric character 
17 1 . In reference 1^ it is shown that the JSD can be taken as a unifying distance between probability distributions. 

In our previous study of the quantum JSD we showed that it verifies all the properties required for a good measure 
of distinguishability between quantum states. In this paper we investigate the metric property of the quantum 
JSD (QJSD), that could be regarded as essential to check on the convergence of iterative algorithms in quantum 
computation [31 . 

The structure of this paper is as follows: next Section is devoted to the formal definition of the classical and QJSD. 
In Section III we investigate the metric character of the QJSD. In the first place we consider the pure states case 
and then we investigate the metric properties for arbitrary mixed states recourse to numerical simulations in different 
Hilbert spaces. Finally, some conclusions are drawn in Sect IV. 

II. CLASSICAL JENSEN-SHANNON DIVERGENCE AND ITS QUANTUM EXTENSION 

The classical JSD between two (discrete) probability distributions P = {pi,p2, ■■■,pn) and Q — ((/i, (72, •■•9w), 
SiPi = 2i ft = 1 is defined as 



(6) 
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where S{P, Q) = J^iPi ^^S Kullback-Leibler divergence. Djs{P, Q) can be also expressed in the form 



1 
2 



P + Q 



log 



-H{P) - -H{Q) 



Pi 



51 * 



2q^ 



Pi + ft 



(7) 



where H{P) = — J^iPi ^ogPi is the Shannon entropy. The classical JSD exhibits several interesting properties. Among 
them we recall the following ones 

• Djs{PtQ) is symmetric and always well defined; 

• it is bounded 

0<Djs{P,Q) < 1, 

and, as it was already stated, 

• its square root. 



djsiP,Q) = y/Djs{P,Q) 
verifies the triangle inequality Eq. ^ (but Djs does not). 



(8) 



A proof of this last fact can be found in references [l7|, ll8| . Alternatively, this can be proved by using some results 
of harmonic analysis due to I. Schoenberg [l^. The basic fact that makes Schoenberg's theorem applicable to the 
classical JSD resides in that it is a definite negative kernel, that is, for all finite collection of real numbers {ci)i<N, 
and for all corresponding probability distributions (-Pi)i<Arj the implication 



N 



J2c^^0^Yl C^CJDJs{P^, Pj) < 



(9) 



i=l 



is valid. A corollary of Schoenberg's theorem allows one to assert that the probability distributions-space, with the 
metric ([U, can be isometrically mapped into a subset of a Hilbert space [20| . 

The classical JSD can be used to distinguish two probability distributions and therefore can be used as well to do 
so for two quantum states described by their density operators, say, p and a. Indeed, let us suppose we choose a 



positive operator value measure (POVM), — that generates two probability distributions via 



P^ 



Tr{E,p) 
rr(E,a), 



for i = 1, ...,M. Then we can use the JSD (O to distinguish between these two distributions. In this procedure we 
have the freedom of choosing the POVM which most clearly distinguishes pi from g;, that is, which makes the value 
of Djs{pi, Qi) the largest. This reasoning motivates to introduce the quantity 



-Djsi(p, cr) = snp Djs{pi,qi 



(10) 



where the supremum is taken over all POVM's. Physically Djsi gives the best discrimination between the states p 
and a that we can achieve by means of measurements. 

By mimicking the extension of Kullback-Leibler divergence to the realm of quantum theory, we define the Q JSD as 



Djs{p,a) = ^ 



S[p 



P + (T 



S a. 



P + (T 



2 / V 2 

that can be recast in terms of the von Neumann entropy Hn{p) = —Tr(p\ogp) in the fashion 

Djsip,a) = Hn (^) - \hn{p) - \hn{o), 



(11) 



(12) 
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This quantity is always well defined, symmetric, positive definite and bounded (0 < Djs{p,<j) < 1). By using the 
corresponding properties of the relative entropy |2l| and expression (jlip it can be shown that, for arbitrary p and tr, 
the following inequality 



Djsip,(j) > Djsi{p,(t), 



(13) 



is valid. The equality is satisfied if and only if p and a commute, that is, the upper bound in (|13p is, in general, not 
attainable for any POVM. 

To conclude this section we give the explicit expression for the QJSD in terms of the eigenvalues and eigenvectors 
of the operators involved in its expression. 



D,js{p,cr) 



1 



^ I (tfck.) prolog 



2r, 

Tk 



Y,\{tk\s,)\^s,\og 



2s, 

Tk 



(14) 



where p = Y.i 'ri\ri){ri\, a = Y,i Si\si){si\, {p + a) ^ I]^»l^i>(*i| and Tk = J2i'^i\{tk\ri)\'^ + I]i SiK^fclsi)!^ 

It should be noted that, when p and a do not commute, the structure of ([M]) is quite different from that of ([7]). 



III. THE METRIC CHARACTER OF THE QUANTUM 

In this section we investigate the putative metric character of the QJSD, that is we try to ascertain whether the 
square root of the QJSD, 

djs{p,<j) - VA/s(p,fT) (15) 

verifies the triangle inequality. The other three properties for a metric are obviously verified by (|15|) . A formal proof 
of property ([3]) for Djs{p, a) has until now eluded us. Unfortunately there is no analog of Schoenberg's theorem 
when operators are involved. Still more, there is no direct way of verifying condition ([9]) for expression (|14p. No 
extension to the case of the QJSD of the proof given in Tt*] has been possible. Incidentally it should be observed that, 
if the upper bound in could be attained for some POVM, the proof of the triangle inequality for would be 
obvious (because y/Djsi verifies it). 

The results to be presented here correspond to a separate analysis of the metric condition for (fT5|) for the two cases: 
when is restricted to pure states and when it acts on the complete set B{H)i . In the first instance we were able 
to give a formal proof of inequality ([3]); in the second one, we checked it by means of a numerical algorithm. 



1. Pure states 



For a pure state the von Neumann entropy vanishes. Then, for two pure states, 

and a^\ip){(f\, (16) 

the QJSD ([12]), becomes 



Djsip,a)^HNl^^j (17) 

After some algebra, we can rewrite (|17|) in terms of the inner product ( | ) : 
Djs{p,<7) = $(|(V'|^)|) = 

The entropy of the average ^ (|'0)('!/'| + can be interpreted to the light of quantum information theory. 

Indeed, let us suppose that Alice has a source of pure qubit signal states lip) and \(p). Each emission is chosen to be 
IV') or \(p) with an equal prior probability a half. Then the density matrix of the source is ^ (|'0)(V'l + Alice 
may communicate the sequence of states to Bob by transmitting one qubit per emitted state. But according to the 
quantum source coding theorem, (|17p gives the lowest number of qubits per states that Alice needs to communicate 
the quantum information (with arbitrarily high fidelity) [23 |. 
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Let us take two fixed arbitrary pure states p = IV") ('01 and u = \<p){'p\ and an arbitrary third one, ^ = |x) (xl- Denote 
the absolute value of the inner products |('0|(p)|, KV'lx)! a-^d l(xl'^)l with x,y and z, respectively, and introduce then 
the function 

In terms of these variables the triangle inequality for (|15|) reads: 

0<G{x,y,z). (19) 

We can decompose the vector \x) into i) a part belonging to the plane determined by and \ip) and ii) another part 
perpendicular to that plane: 

\x)=a\iP)+b\ip) + \xi_). 

with \a\ < 1 and |6| < 1. Then 

y = \a + b{ilj\ip)\ and z ^ \a*{(p\ip) +b*\. 

As a function of a and b, for x fixed, G is a concave function on the circles \a\ < 1 and \b\ < 1 (in the sense that its 
second derivative is negative) and it vanishes for y = x and z = x. This guarantees that inequality (|19p is satisfied 
for arbitrary y and z. 

2. Arbitrary mixed states 

Here we attempt a numerical verification of the triangle inequality for the distance (|15p when arbitrary mixed 
states are involved. As a first approach, we numerically evaluate the inequality ([3]) by generating random states in a 
N-dimensional Hilbert space. The space of all (pure and mixed) such states can be regarded as a product space of 
the form [H: 

n^v X A, 

where V stands for the family of all complete sets of ortho-normal projectors {Pi}f^ , J2i A = I (I the identity matrix), 
and A is the set of all real A^— tuples of the form {Ai, . . . , Xn}', K G R; J^i Xi ^ 1; < < 1. Any state in H is of 
the form p = A^Pi. 

In exploring exhaustively Ti. we need to introduce an appropriate measure p on this space. Such a measure is 
required to compute volumes within 7i, as well as to determine what is to be understood by a uniform distribution of 
states on Ti.. The measure that we adopt here is taken from the work of Zyczkowski et al. (25., ,26]. 

An arbitrary (pure or mixed) state p of a quantum system described by an A^-dimensional Hilbert space can always 
be expressed as a product of the form: 

p = UD[{\i}]U'^. (20) 

Here U \saivN xN unitary matrix and Z)[{ A^}] is an A^ x A^ diagonal matrix whose diagonal elements are, precisely, our 
above defined {Ai, . . . , \n\- The group of unitary matrices U{N) is endowed with a unique, uniform measure, known 
as the Haar's measure, v [23|. On the other hand, the A^-simplex A, consisting of all the real A^-uples {Ai, . . . , Ajv} 
appearing in (PO)) . is a subset of a (A^ — l)-dimensional hyperplane of R^. Consequently, the standard normalized 
Lebesgue measure Cn-i on R^~^ provides a measure for A. The aforementioned measures on U{N) and A lead to 
a measure p on the set S of all the states of our quantum system [2^, |2^, [23| , namely, 

fi = u X Cn-1- (21) 

In our numerical computations we randomly generate mixed states according to the measure (|2T|) . In order to assess, 
for these randomly generated states, how the triangle inequality ([3]) is satisfied, we define the auxiliary quantity 



Adjs(p, cr) = djs{p, + djs{^, cr) - djsip, <j) 



(22) 
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and evaluate it for a large enough number of simulated states. This procedure is repeated for different dimensions of 
the Hilbert space. 

We investigate the positivity of Adjs, upon which the metric character of the square root of the QJSD is based, 
by constructing the probability distributions V for the values of Adj^. The corresponding histograms, for different 
dimensions of the Hilbert space, are depicted in Fig. 1. As we are mainly interested in the positivity of Adjs, we 
just plot the tails of the concomitant distributions V, selecting the portion for which one has, say, Adjs < 0.2. Such 
a choice allows us to portray in sufficient detail the region of the distribution where a violation of the inequality ^ 
can be detected. 

The probability for the particular value AcJjs = actually represents the probability for finding a triplet of density 
matrices for which Adjs < 0. None such triplet of states has been found, which entails that the probabilities 
for violating the triangular inequality vanish for all the distinct Hilbert-space dimensions we have considered here. 
Actually, the probability for low values of Adjg becomes significantly smaller as the dimension of the pertinent Hilbert 
space under study augments (the PDF's for higher Hilbert space dimensions than those here reported have been also 
computed). 

The total number of randomly generated states was rather large (10*) in order to obtain a sufficiently large number 
of points belonging to the tail-regions. These points fall then within the zone of low probabilities. The fact that 
no triplet of states violating inequality ([3]) has been encountered could be thought of as being numerical evidence 
for the metric character of the square root of the QJSD. The distributions in Fig. 1 clearly depend on the measure 
(PT|) used to compute them. Higher probabilities for low values of Adjs can actually be obtained if one restricts 
the computation of the histograms to states with a high degree of mixedness, although it must be noted that such 
probabilities still diminish as the dimension of the associated Hilbert space grows. 

To avoid a statistical dependence on the measure (|2ip we propose an alternative numerical approach by performing 
a numerical minimization of Adjs- Any quantum mixed state is completely determined by a finite number Hp of 
parameters which depends on the dimension of the Hilbert space. To determine the minimum possible value of Adjs , 
one needs to find the optimal values for such parameters. To such an end we use a simulated annealing algorithm in 
which the parameters are iteratively modified until convergence to the optimal values is reached. 

After running this algorithm for different Hilbert space dimensions and for different triplets of initials states, one 
detects always convergence to the same solution 



mm{Adjs{p,t'j)}^0. (23) 

The optimal situation is reached when p and ^ are equal. In our numerical search these states are always found to 
coincide with the maximally mixed state for the Hilbert space-dimension considered in each case. It is actually not 
enough to minimize Eq. (|22p . because we wish it to be a minimum for any of the three different ways to order the 
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three states. If we minimize the average of those three possible orderings, the minimum is also Adjs{p, <j) = 0, and 
it is obtained when the three states become the maximally mixed state. 

This last method, although does not provide us with a formal proof of the metric character of the square root of 
the QJSD for mixed states, does yield a clear and strong evidence about the validity of the conjecture advanced in 
the initial part of this paper, that constitutes the leitmotif of this work. 



IV. CONCLUSIONS 



The main purpose of this work was to investigate the metrical property of the QJSD. We were able to show that the 
square root of the QJSD verifies the triangle inequality, giving to this distance the character of a metric. Although we 
have proved this claim (for mixed states) only by giving numerical evidences, we believe that the cases here analyzed 
are sufficiently representative so as to render credible the claim that metric properties are verified in general for the 
QJSD. 

A second item deserves to be pointed out, which emerges from the following two facts: 

• On the one hand, we have showed that, when restricted to pure states, the square root of the entropy of the 
average ^ {\ip){ip\ + \ip){Lp\) is a true metric. 

• On the other hand, a classical result from Uhlmann Q asserts that the fidelity of states p and a 



Fip,cr) ^ Tr\J p^ap-2 

can be expressed in the form 

max KV'lv')! (24) 

where the maximization is over all purifications j?/;) of p and all purifications \f) of ct [28| . 

These two facts motivate us to introduce an alternative metric for arbitrary mixed states. Given two arbitrary 
mixed states p and a we can define 




dnip, cr) = , mill \ Hn [ (25) 



where the minimum is taken over all purification j?/;) of p and all purifications \ip) of a. In (|25p we must look for the 
minimum, not for the maximum as in (j24p . due to the decreasing nature of $, eq. (|18p . as a function of 

Obviously the basic properties required for a good distinguishability measure are inherited by from those 
verified by the QJSD. Additionally, several interesting questions arise from this proposal. For example, what relations 
exist between (l25t and (fTTj) : or, in general, how ((25)) relates to other quantum distances. A more detailed study of 
the properties of this quantity will be presented elsewhere. 
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